The magnetic susceptibility of disordered non-diffusive mesoscopic systems 
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Disorder- induced spectral correlations of mesoscopic quantum systems in the non-diffusive regime 
and their effect on the magnetic susceptibility are studied. We perform impurity averaging for non- 
translational invariant systems by combining a diagrammatic perturbative approach with semiclas- 
sical techniques. This allows us to study the entire range from clean to diffusive systems. As an 
application we consider the magnetic response of non-interacting electrons in microstructures in the 
presence of weak disorder. We show that in the ballistic case (elastic mean free path £ larger than 
the system size) there exist two distinct regimes of behaviour depending on the relative magnitudes 
of I and an inelastic scattering length L^. We present numerical results for square billiards and 
derive approximate analytical results for generic chaotic geometries. The magnetic field dependence 
and dependence of the disorder-induced susceptibility is qualitatively similar in both types of 
geometry. 
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I. INTRODUCTION 

Mesoscopic physics has traditionally involved the 
study of nwel phenomena in phase-coherent, disordered 
conductoril. In the diffusive regime, where the elastic 
mean free path £ is much smaller than the system size 
L, the electron motion resembles that of a random walk 
between impurities. 

The more recent development of high-mobility semi- 
conductor heterostructures together with advanced litho- 
graphic techniques has led to the confinement of elec- 
trons in two-dimensional microstructures of controllable, 
non-random geometry. Thereby, a further regime, coined 
"ballistic" since £ > L has been realised. Such technolog- 
ical achievements have motivated theoretical approaches 
where the actual microcavities are approximated by 
"clean" quantum billiards ignoring impurity scattering 
completely. In these models the electron motion is only 
affected by bounces at the confinement potential, which 
is the opposite case to the diffusive regime where con- 
finement effects are not important on timescales shorter 
than the Thouless time. 

However residual impurity scattering is nearly un- 
avoidable even in high-mobility microstructures. It has 
been shown that weak disorder can be strong -ejiiough to 
mix energy levels, influence spectral statisticsBu and af- 
fect related thermodynamic quantitiesQ in the "ballistic" 
regime. 

In defining spectral correlation functions in the "bal- 
listic" regime one has to distinguish between disorder av- 
eraging, which we denote by (. . .)^, and size averaging, 
(. . .)^, which we assume to be equivalent to energy av- 
eraging. Pure disorder averaging corresponds to the ex- 
perimental situation of an ensemble of weakly disordered 
microstructures of the same size. Recent work has shown 
that for this case spectral correlation functions contain, 
and often are dominated by, strong oscillatory structures 
in energy (on scales kpL, kp being the Fermi wavenum- 



ber) reflecting the presence of the confincmcntao. They 
can be semiclassically interpreted as density-of-states os- 
cillations related to classical periodic orbits of the corre- 
sponding clean system. 

In the present work we consider spectral correlations 
after both energy and (independent) disorder averaging 
corresponding to the experimental situation of an ensem- 
ble of disordered systems with variation in their sizes. 
Then one is able to divide the two-level oorrelation func- 
tion, K{ei,e2)^ into two separate terms,B'El 

Kiei,S2;H) = {K\si,e2;H))^ + K''{suS2;H). (1) 



Here 

K''{si,e2) 
K''{ei,e2) 



Meiy{e2))d-He,)),He2))Alv\ (2a) 
{Hei))dHe2))d)L/y^-l, (2b) 



where v denotes the single particle density of states 
and V = {{i'{e))d)L its mean part. K"^ is a measure 
of disorder-induced correlations of the density of states, 
while is given by size-induced correlations. In a dif- 
fusive system, £ < L, the disorder-averaged density of 
states {v{e)) J is a constant, so that is vanishingly 
small and K dominates. However, for £ > L the density 
of states contains terms which oscillate like cos{kpL), and 
both correlation functions may be relevant. In contrast, 
once £ > L{kFL)'^~^ disorder-induced mixing of levels is 
negligible and ^ prevails. 

The orbital magnetism of mesoscopic quantum sys- 
tems is sensitive to spectral correlations as e.g. mea- 
sured by K{ei,S2) and therefore has been the subject 
of much theoretical interest^ as well as experimental 
investigationsOiij. For isolated systems with a flxed 
number of particles it is necepaHt to consider averag- 
ing under canonical conditionaHTLJ resulting in a large 
contribution to the average magnetic response which can 
be by orders of magnitude larger than the (bulk) Landau 



1 



diamagnetism. The corresponding susceptibihty is given 
by0 



2 j dm 



(3) 



for temperatures k^T larger than the mean level spacing 
A. In Eq. (|) H is the magnetic field and {6N^{^i;H)) 
is the variance in the number of energy levels within an 
energy interval of width equal to the chemical potential, 
fi. This variance is related to K{ei,e2]H), Eq. by 
integration of the level energies £i, £2 over the energy 
interval. Throughout this work we will label the con- 
tributions to the susceptibility, corresponding to ( K"^ ) ^ 
and K^, as (x''(-ff)) and respectively. Hence 

the total orbital magnetic susceptibility is composed of 



(4) 



Diagrammatic techniques to treat impurity scatter- 
ing are usually designed for diffusive systems, £ < L, 
and therefore rely on the assumption of translational 
invariance. In this paper we combine a diagrammatic 
perturbation approach with semiclassical techniques to 
perform energy and impurity averaging and to calcu- 
late the disorder induced part of the energy correlation 
function and the related susceptibility (x'^(7J)) for disor- 
dered microstructures in the non-diffusive regime. In this 
quantum-semiclassical hybrid approach, scattering at the 
impurity potentials, which are assumed to be (5-like, is 
treated quantum mechanically in a perturbation series 
with respect to the disorder. Boundary effects are incor- 
porated in a semiclassical representation of the Green 
functions, which enter into the impurity diagrams, in 
terms of classical paths. This method, therefore, takes 
into account, in a systematic way, contributions from 
(closed) trajectories which involve both scattering at im- 
purities and specular reflection at the confinement po- 
tential. The procedure allows us to study the complete 
crossover from diffusive to clean systems. A preliminary 
brief account of this work has been presented in Ref.El. 

Our approach carries features of the "method of tra- 
jectojies" originally devised for thin supepeonducting 
filmstil. Note that similarEJ, and alternativeEj, methods 
have been also employed to consider weak localization in 
thin films in a parallel magnetic field. Related semiclas- 
sical approaches for diffusive systems were proposed by 
Chakravarty and Schmidtfl aiid Argaman et al.Ej. More 
recently, Agam and FishmanEZI studied the spectral form 
factor for ballistic systems with rigid disks (spheres) as 
impurities from a quantum chaos point of view. We fur- 
ther note that a non-perturbative approach for ballistic 
systems has been developed by; Muzykantskii and Khmel- 
nitskii, the ballistic a modelE3. 

For real systems, besides £, there are additional rel- 
evant lengthscales at which inelastic scattering (i^) or 
temperature smearing (Lt) produce a damping of prop- 
agation. For clarity we refer to such a lengthscale in the 



following as L^, although we assume that similar general 
arguments will hold for finite Lt- For ballistic motion 
at the Fermi energy Ep, we can relate to the level 
broadening 7 by 



T 



kpL A 
27r 7 



(5) 



It divides the "ballistic" regime into two sub-regimes. In 
the first, L, < £, the particle motion is nearly ballistic: 
the damping (of the Green function) due to inelastic scat- 
tering typically occurs before elastic impurity scattering; 
for the remainder of this paper we refer to this regime as 
inelastic. In the second regime, L < £ < L^, a particle 
may scatter many times off impurities before scattering 
inelastically. We call this regime elastic. In our semiclas- 
sical treatment we consider the contribution of orbits of 
all lengths (smaller than vpt^, where tu ~ fi/IS. is the 
Heisenberg time) in both the elastic and inelastic regime. 

We present and compare results for non-interacting 
ballistic quantum systems with integrable and chaotic 
classical dynamics in the clean limit. We derive ana- 
lytical estimates for the average susceptibility (x'^(iJ)) 
of generic chaotic microstructures assuming ergodicity 
for the classical paths involved. As an example of an 
integrable geometry, we treat in detail the repiresenta- 
tive case of the square billiard. Experimentsllj on the 
magnetic susceptibility of ensembles of squares were per- 
formed in the "ballistic" inelastic regime, mo1jij\pi;ing the- 
oretical studies of the susceptibililaiifor L < ffl'EJ. Gefen, 
Braun, and Montambaux (GBM)cil considered the con- 
tribution of trajectories longer than £ to (x'*(^)) in 
approximate way, finding a paramagnetic ^-independent 
contribution at|-pero field, whereas Richter, UUmo, and 
Jalabert (RUJ)^ calculated {x^{H)) for a square by as- 
suming that the disorder perturbs the phase, but not 
the trajectory, of semiclassical paths of the corresponding 
clean geometry. In the elastic regime we find agreement 
with the results of GBM with regard to the kpL behav- 
ior, while the zero field susceptibility is weakly dependent 
on I. In the inelastic regime, for larger level broadening, 
we find entirely different results, namely an exponential 
dependence on L/L^ and on £/ L. 

The present paper is organized as follows: In the next 
section we outline our combined semiclassical diagram- 
matic technique. In Sec. Ill we illustrate our approach 
by applying it to the case of integrable (square) billiards. 
We give analytical results for spectral correlation func- 
tions (at zero magnetic field) and numerical results for 
the magnetic susceptibility. In Sec. IV we present the 
derivation of the corresponding disorder induced spec- 
tral correlations and susceptibility for chaotic ballistic 
geometries. 
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II. SEMICLASSICAL DIAGRAMMATIC 
APPROACH 

In this section we summarize our semiclassical evalu- 



ation of the disorder correlation function K'^^ Eq. (2a), 
and the corresponding susceptibility (x^iH))- To this 
end we begin with the diagrammatic formulation of the 
problem and then perform the semiclassical approxima- 
tions. 



schematically in Fig. ||a. Note that the sum of diagrams 
sl^^ is equivalent to a number of one-loop diagrams in 
the conventional notation: they include the dominant 
contiibution to K''^ in the diffusive regime for d > 2 
(RefH3) but are actually smaller than some two-loop dia- 
grams in a wide region of energies in the diffusive regime 
at d = 2 (Rcf .E^) . However in the region of interest of this 
paper, the non-diffusive "ballistic" regime, the diagrams 
sIt^^ dominate. 



A. Diagrammatic framework 

We consider non-interacting electrons in a weak, per- 
pendicular magnetic field and a random Gaussian poten- 
tial V , with {V{y)) — and correlator 



{V{v)V{v'))^ 



6{r - r') , 



(6) 



describing white noise disorder. In Eq. r is the mean 
elastic scattering time, £ = vpT, and D = mL^ /2Tih^ 
in two dimensions. In terms of the retarded and 
advanced single particle Green functions obeying the 
boundary conditions of the corresponding clean system, 
r2\e \ H), the correlator K'^ may be written as 

(7) 

Here, the average is taken over impurities only (for a 
given system size) and the symbol ((. . .))^ implies the in- 
clusion of connected diagrams only. 

We are particularly interested in the field sensitive 
part, K"^, of A'''.!-. Following a diagrammatic approach 
introduced in Rcf.H it can be expressed as 



(C) 

with uj = El — £2- The Cooperon type diagrams Sn are 
defined by 

5(^)(c.;i/) = Tr[c(^)(c^;iJ)]" (9) 
/ n«^''-^ n C^''Hr™,r™+i;c.;iJ) 

/ i=l / rn=l 

with r„+i = ri and 

C(^)(ri,r2;cc>;i/) = 

-^G+(ri,r2;£i;i/) G- {r,,r2;e2; H) . (10) 
Zni'T 

Here G^^^^ = {G^^^^)^ is the disorder averaged single 

(C) 

particle Green function. An example, 54 , is shown 




(a) 



(b) 



FIG. |l| The diagrams si^\ (a) shows a schematic form of 



the diagram 51 , (b) shows a pair of typical real space semi- 

(c) 

classical trajectories which contribute to 54 . 



B. Semiclassical treatment 

Due to the lack of translational invariance in confined 
systems it is no longer convenient to evaluate diagrams 
such as that in Fig. ^ in momentum space. Instead we 
work in configuration space and compute the integrals 
(^) invoking a semiclassical approximation. 

Semiclassically, the single particle Green function 
Q'^(ri,r2) can be expressed jas a sum over classical tra- 
jectories t between ri and r2Ea. After disorder average it 
readsQ 

G+(ri,r2)^ J2 Acxpfj5,(ri,r2)-§) . (11) 



t:ri-*r2 



The classical amplitude Dt includes the local density of 
trajectories near the path St stands for the classical 
action along the orbit including the Maslov index, and 
Lt is the orbit length. 

Eq. (|ll] ) was derived in a semiclassical framework un- 
der the assumption that weak disorder modifies merely 
the phases St, which leads to damping on the scale of £, 
leaving the trajectories unaffected; i.e. the sum is taken 
over the paths of the corresponding clean system. 

For the bulk case this treatment is equiivalent to the 
usual diagrammatic treatment of disorderQ: For white 
noise the self-consistent Born approximation leads to the 
following integral equation ior the impurity averaged sin- 
gle particle Green functioncJ: 

G+(ri,r2) =^+(ri,r2)+ (12) 
+ ^ J dra ^+(ri,r3)G+(r3,r3)G+(r3,r2) . 



3 



In the bulk case G+(r3, r^) in the above equation is given 
in terms of "paths of zero length" leading to the expo- 
nential damping of the (free) Green function without dis- 
order. For confined systems further paths of finite length 
may semiclassically contribute to G+(r3,r3). They can 
be viewed as paths scattered off the impurities which 
may lead to additional corrections to the impurity av- 
eraged single particle Green function. Howejier, these 
contributions are of higher order in h and 1 /tvB. 

Upon using the semiclassical expression (^l|) in 
Eq. (|lO|) the two-particle operator (^^\ri, uj; H) is 
then given in terms of a double sum over pairs of classi- 
cal paths which explicitly include the effect of boundary 
scattering. However most pairs (of different paths) pro- 
duce oscillating contributions-which we assume to vanish 
after energy or size averagingEH. Therefore the main con- 
tribution to the field sensitive part of ( K"^ )^ arises from 
diagonal terms (otherwise known as the Cooperon chan- 
nel) obtained by pairing paths with their time reverse. 
Assuming that the magnetic field affects the phase of the 
particles but not their trajectories we can write for the 
actions 



^5t(e,;i/)~ 

lSt{EF;H^Q) + {e,-EF)Tt + — / ' A.dr , (13) 

where Tt is the period of the trajectory, A is the vector 
potential and ipQ = hc/e. The diagonal approximation 
for C'*^-* then gives 



C(^)(ri,r2;c.;ff) 



Cf'(ri,r2;^;F) , (14) 



where 



X exp 



Lt 



VF\Dt\' 

-f iu)Tt 4- 



.47r 



A.dr 



(15) 



The level broadening (implicit in i^, Eq. (||)) was in- 
troduced via uj ^ uj + i^. Eq. (^ depends, besides 
only on the system without disorder and holds for both 
integrable and chaotic geometries. 



C. Disorder induced magnetic response 

Using Eqs. (^) and (^ the disorder induced contribu- 
tion to the average magnetic susceptibility is given by 
((^ = i/LV^o) 



\XL\ 



^2 Q^2 Z_-/ ji 

71 — 1 



^V)(0;^), (16) 



where the bulk Landau susceptibility is xl 
— e^/247rmc^ for spinless electrons. 



Assuming that £ remains fixed as kpL changes, one 
sees from Eqs. (jl|) and (|l|) that {x'^[ip)) contains the 
dimensionless variables 7/ A and kpL only in the com- 
bination / l^)/kpL, since they are absorbed into the 
dimensionless variable L^/ Eq. (|^). It is now assumed 

that the in Eq. ( |l6| ) contain diagonal terms only. 

Therefore the propagator C^'-^^ (Eq. (p^) is made up of 
a summation over all diagonal pairs of paths including 
boundary scattering between any two given impurities 
situated at ri and V2. On taking the trace over n prop- 
agators C'"'~^\ one sees that the field sensitive part of S'„, 
Eq. (H) , consists of a summation over flux-cnclosing pairs 
of closed paths (in position space) involving n impurities 
and an arbitrary number of boundary scattering events. 
An example of a pair of paths contributing to ^4 is 
shown in Fig. |^. 



III. INTEGRABLE GEOMETRIES: 
SQUARE BILLIARDS 

In the following we apply the above formalism to com- 
pute the magnetic response of an ensemble of disordered 
billiards with regular geometry. We illustrate the method 
and present numerical results for the case of square bil- 
liards. 



A. Numerical technique 

We consider specular reflection at the ^-boundaries 
and employ the extended zone schemaHd'E^ to write 
C,^^\v I, V2\0J -,11), (Eq. (p^), as a sum of propagators 

along straight line pathsQ^ (vi^v\\uj\H)^ of the form 
Eq. (p^), where v\ are images of the position V2- The 

diagrams 5^*"^ (Eq. (||)) are then calculated by diagonal- 
ising C,'^'-^\ 

At zero magnetic field this approach enables one to re- 
cover the result of Ref .a for the approximate diagonalisa- 
tion of the spectral correlation function K'^iei , 62; H = 0) 
as outlined in the Appendix. 

For finite magnetic field the integrals over the mag- 
netic vector potential along the paths do not allow for 
an analytical diagonalisation of ^'^'-^K However we are 
able to use the fact that all the variations of C'"'^^ occur 
on classical lengthscales; rapid oscillations on the scale 
of Xp cancel out. This enables an efficient numerical 
computation. To this end we discretise the configuration 
space of the square billiard using a lattice with grid size 
greater than Xp- By summing over all trajectories (up to 
a length ^ L^) which connect two lattice cells we com- 
pute the corresponding matrix elements of C^*^^ in this 
representatiorL_Aftcr diagonalization we obtain (x'^(i?)) 
from Eq. (p^EEl. This method is not restricted to the 
square geometry but can be in principle applied to any 
geometry. 



4 



(a) 






FIG. y. Examples of diagonal pairs of real-space trajecto- 
ries with one impurity scattering position in a square, (a), 
(c), (d) show pairs of typical recurrent orbits with momen- 
tum exchange at the impurity which contribute to <^x'^(^))i 
whereas (b) shows a pair of orbits following periodic orbits of 
the corresponding clean system which contribute to <^x^(^^))- 

(C) 

Note that the diagrams Sn do not include closed pairs 
of paths which follow periodic orbits of the correspond- 
ing clean system. Such paths involve zero momentum 
transfer between the Green functions at the impurity po- 
sitions. They actually represent disconnected diagrams 
and are included in ^x^(iJ)^ which has been computed 
in Ref.Q. It is the presence of these periodic orbits in the 
determination of {x^{H)'^ that leads to strong sensitiv- 
ity of the-|(Mhi-tal magnetism with respect to the system 
geometryQ'CJ'Lll. The above numerical procedure does not 
distinguish between connected contributions to 5i, with 
momentum exchange at the impurity as in Figs. ^,c,d, 
and disconnected diagrams involving one common impu- 
rity position. The dominant disconnected diagrams arise 
from the shortest flux enclosing periodic orbits of length 
Lt — 2^/2L, Fig. and their repetitions. An estimate 
of thfip,magnitude is obtained following the approach of 
RUJltEj. On summation over all repetitions of the fun- 
damental orbit we find 



(x(o)> 



dis 



1 



\Xl\ 



8V2L 

5tt £ smh^[V2{L/e + L/L^)] 



(17) 



To avoid double counting we explicitly subtract this es- 
timation from the numerical determination of ^x'^(-f^))- 



B. Zero field susceptibility 

We briefly summarize our results for the magnetiq-sus- 
ceptibility at zero field and refer the reader to Ref.B for 



further details. The technique covers the whole range 
from the diffusive regime to the clean limit and yields an 
averaged susceptibility which is always paramagnetic. In 
the diffusive limit, £ < L, the disorder induced suscepti- 
bility increases linearly with £ in agreement with Ref.E3. 
In the ballistic regime, (x'^(O)) exhibits a maximum as 
a function oi £/ L. The occurence of the maximum may 
be related to the competition between different effects of 
the impurity scattering on C,^^^ , Eq. (|l5|): while the single 
particle Green functions are exponentially damped with 
£, it enters as £^^ into the prefactor. 




FIG. ^ Comparison of the relative contributions of diagrams 
Si"^' to (x''(0)) for l/L = 4 as a function of L^/L. The solid 
line corresponds to the inclusion of all diagrams while, from 
the bottom, the dotted line corresponds to only, the 

dashed-dotted line to S{ 



line to S^^^ and 53'"'''. The long dashed line is the con- 

tribution of (x^(0)). 

In the elastic regime, L < £ < L^, (x'^(O)) ex- 
hibits a weak dependence on £ hut is on the whole 
close to the prediction by GBMcJ, which is a para- 
magnetic ^-independent contribution (x'*(0)) / |xl| ~ 
0.23 kFL{A/'j), For larger £, in the inelastic regime, 
L, < £, the susceptibility decays exponentially with 
both £/L and L/L^. 

The disorder induced susceptibility {x'^{0)) for £/L = 
4 and a range of /L including both the elastic and in- 
elastic regimes is shown as the solid line in Fig. We also 
include the contributions of diagrams Eq. (^), with 

a different number n of impurity scattering events in or- 
der to analyse their relative weights. From the bottom, 

(c) 

the dotted line corresponds to only, the dashed- 



;(C) 



(C) 

and 52 , and the short dashed 



dotted line to Sl*^^ and 82^' , and the short dashed line to 



s['"\ 82^^ and S^^' . For completeness, the long dashed 
line is the contribution of (x'^(O)) (see below). The solid 
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line curve shows linear behaviour in the elastic regime, 
as predicted by GBM, and exponential behaviour in the 
inelastic regime. In the limit <^ i, where (x''(0)) 
is very small, the lowest n diagrams dominate but for 
> £ a number of diagrams contribute significantly. 
We note that, contrary to what is stated in Reffl, the 

(C) 

contribution from 5J is not particularly small through 
the whole range of L/L^. 



C. Finite field susceptibility 



We compute the finite field susceptibility for 
weak fields, assuming that the field does not perturb tra- 
jectories away from their zero- field straight line paths, i.e. 
we assume rc > i where Vc = Tikpc/eH is the cyclotron 
radius. Fig. || shows (^x'^{H)'j for kpL — 60 and 7/A = 1 
{L^/L — 9.55). From the top (at H — 0) the mean free 
path is ^/L = 5, 10, 20 and 50. As the field increases from 
H = the susceptibility falls in magnitude and eventu- 
ally becomes diamagnetic. The field value at which the 
susceptibility crosses over to diamagnetic behaviour de- 
creases as i/L increases. This implies that the typical 
area enclosed by a pair of recurrent orbits increases with 
i. 

In the range of validity of the above approach, £ < 
Tc, kpL/^ and 7 A we do not observe oscillations of 
(x'^{H)) as a function of field as seen by GBMcil (albeit 
for a smaller value of 7/A). 



2^ 
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FIG. 0. Field dependent susceptibility {x'^iH)'^ for a square 
for kpL = 60 and 7/A = 1 {L^/L = 9.55). From the top (at 
H = 0) the mean free path is £/L = 5, 10, 20 and 50. 



D. Comparison of the combined semiclassical 
contributions with quantum mechanical results 



In this section we compare semiclassical results for 
(x'^(-ff)) and (x^(^f))|-|With numerical quantum calcu- 
lations taken from Ref.u. In order to obtain an analytic 
expression for|-|the contribution of size-induC|ei-correla- 
tions (x^(iJ))la, we employ the results of RUJa-Ea at zero 
temperature and introduce the level broadening 7 in the 
same wjav as for the disorder correlations. It has been 
shownEli23 that the low field susceptibility of an ensem- 
ble of clean squares is dominated by the shortest flux 
enclosing periodic orbits of length Lt = 2\/2L and their 
repetitions over a broad range of temperature (or corre- 
spondingly inelastic scattering strengths). For the ballis- 
tic white noise case considered here, the effect of disorder 
averaging on the susceptibility was described by an addi- 
tional damping exp{—Lt/i) of the response of the clean 
systemQ. This result corresponds to (x^(^)) including 
the disorder damping exp(— Lt/2£) of the two single par- 
ticle Green functions. 




FIG. 



i 



Disorder- and size-induced contributions to 



the susceptibility. The mean susceptibility for small field, 
ip/ipo — 0.15, and high temperature fcsT/A = 2 {Lt/L ~ 1.8 
at kpL = 70), for various disorder strengths, taking into 
account the variation of the mean free path with kp; I oc kp. 
From the top, the elastic mean free path is l/L = 8,4,2,1 
(at kpL « 70). The full and dashed curves are data from 
a quantum mechanical calculation and from a semiclassical 
evaluation of the contribution of size-induced correlations, re- 
spectively, taken from Ref.Q. The dotted curve represents the 
sum of our semiclassical evaluation of the disorder-induced 
correlations, using 7/A « nksT/A, and^the semiclassical 
data for size-induced correlations from Ref.u. 

The numerical quantum calculationsS were obtained 
by diagonalisation of the Hamiltonian for non-interacting 
spinless particles in a square billiard with perpendicu- 
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lar magnetic field and white noise random potential. In 
this calculation temperature smearing, rather than level 
broadening due to inelastic effects, was introduced. 

The full lines in Fig. |^ show the numerical results for 
various disorder strengths as a function of kpL. From 
the top, the elastic mean free path is £/L — 8, 4, 2, 1 (at 
kpL « 70). As kpL changes, the mean free path is as- 
sumed to change also, according to £ oc kp. p. 

The dashed curves in Fig. || are data taken from Ref.cl 
for the semiclassical evaluation of the contribution of size- 
induced correlations to the susceptibility, (x^(i/)). We 
calculated the contribution of the disorder-induced cor- 
relations, (x'*(ff)), using the approximation^ 7/A « 
TrksT/A and taking into account the variation of mean 
free path with kpL. The resulting semiclassical evalua- 
tion of the total susceptibility, (^x^{H)) + {x'^{H)), is ob- 
tained by adding our results to those given by the dashed 
curves; the results are shown as dotted curves in Fig. |^. 

For kpL > 60 the agreement of the semiclassical (dot- 
ted) curves with the quantum results is considerable; it is 
clear that the difference between the semiclassical results 
for {x'^{H)) and the quantum results is accounted for by 
the addition of {x'^{H)). However even for kpL > 60 the 
agreement is not exact; we believe that this minor dis- 
crepancy is due to the different mechanisms of damping 
that were used. 

For kpL ^ 60 there is worse agreement between the 
semiclassical (dotted) curves and the quantum results. 
This may be related to the fact that the semiclassical re- 
sults show no sign of the oscillatory structure present in 
the quantum curves. This structure may be caused by 
off-diagonal pairs of different families of short periodic or- 
bits contributing to (x^(^))i which are not completely 
suppressed through energy average and are not included 
in the present semiclassical approach. 



E. Relation to experiment 

Measurements of the orbital magnetism of ballistic sys- 
tems, which are experimentally, realized as semiconductor 
micipstructures, are still rareQ'EEI. In the experiment of 
Rcf .EiJ the susceptibility of ensembles of squares was stud- 
ied. The samples used had estimated values for the elas- 
tic mean free path oi i/L ~ 1 — 2, for the phase-coherence 
length of (3 — 10)L and for the thermal cutoff length 
of Lt/L ^ 2. Therefore, the Icngthscale (Eq. (^)) is 
determined by the shorter length Lt- 

Fig. H shows that for these experimental parameters 
both the disorder and size-induced correlations are rele- 
vant, however the latter contribution is dominant. 

The above remarks hold for white noise disorder. How^, 
ever, experimental ballistic structures as those of Ref.tHi 
are usually characterised by smooth disorder potentials. 
Smooth disorder effects can be in principle incorporated 
into the present calculation by introducing an angle- 
dependent cross section for the impurity scattering be- 



tween two successive trajectory segments. The effect of 
smooth disorder on (x^) has been analysed in RefJj which 
showed that the reduction of the clean contribution is not 
as strong as that for white noise disorder and no longer 
exponential. We therefore expect that in the parameter 
regime of the experiment the domination of the suscep- 
tibility by size-induced correlations is further enhanced 
when considering smoothed disorder. 

The measured value of the susceptibility at low temper- 
ature was x(0) ~ 100|xl|, with an uncertainty of about 
a factor of four. The combined contributions (x'') and 
(x^), together with an interaction contribution of the 
same ordeiu, are in broad agreement with the experi- 
mental result. We note, however, that a theoretical ex- 
planation of the temperature dependence of the measured 
susceptibility is still lacking. 



IV. CHAOTIC GEOMETRIES 

For systems with a generic chaotic, clean counterpart 
we obtain an analytical estimate for (x''(0)) under cer- 
tain statistical assumptions with respect to the classical 
trajectcpies involved. To this end we make use of the 
relationlij (L^ denotes the system area) 



J2 \D,\H{t-t,)^P{r,,r,;t) 



(18) 



j:ri^r2 



in order to transform the sums over classical densities 
|£>tp in Eq. (|l5|) into classical probabilities P(ri,r2;t) 
to propagate classically between impurities at ri and r2 
in time t. 

In the following we assume that in the "ballis- 
tic" regime £,L^ ^ L the conditional probability 
P(ri, r2; i|^) to accumulate an "area" A during the prop- 
agation from ri and r2 is independent of ri and r2. Fol- 
lowing Ref.Q which considered clean billiards, we find 
from Eqs. ^ and (|l|) 



C(^)(ri,r2;c^;77)~ - / dt d^P(ri,r2; 

Jo J -00 



t\A) 



fATrAH\ 

X cos exp 

\ "Po J 



-jt h icijt 

T 



(19) 



We assume that after several bounces off the bound- 
ary the area distribution P(ri, r2; i|-4) becomes Gaussian 
with a variance a independent of ri and r2 and indepen- 
dent of £. Substituting the above expression for C(^)" into 
Eq. IM) and performing the r, t and A integrals we obtain 



the following closed approximate form for Sn 



5 1- 1 -I- 7T — ILUT 



(20) 



Summation of the 51*"^ Eq. (p!6|), leads to the disorder 
induced contribution to the average susceptibility. 
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\Xl\ 



24 



In 



(21) 



where (p = HL'^/tpo. It is possible to obtain an explicit 
expression for (^x'^i'P)} analytically. For brevity we give 
only the limits. 

The susceptibility has a paramagnetic maximum for 
zero field, 



\xl\ ^L^{L^ + e)' 



(22) 



and it changes sign only once, for example in the elastic 
regime at the critical flux 



L3 



8tt^<j{L^/L) ' 



Ls>e 



(23) 



For large field, ip ipc, the susceptibility remains dia- 
magnetic and decays rapidly. 



6L^ 



\Xl\ 



(24) 
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FIG. ^ Field dependent susceptibility (x''(^)) (normalized) 
for a generic chaotic geometry for kpL = 60 and 7/A = 1 
(L^jL = 9.55). From the top (at H = 0) the mean free path 
is £/L = 5, 10, 20 and 50. 

The susceptibility as a function of field is plotted in 
Fig. I for /cF^ = 60 and 7/A = 1 {L^/L = 9.55). Various 
disorder strengths are shown, from the top (at H — 0) 
£/L — 5, 10, 20 and 50. These curves are qualitatively 
similar to those for the square geometry shown in Fig. ^. 
This may be related to the fact that for trajectories be- 
ing (multiply) scattered at impurities the character of 
the clean geometry, namely regular or chaotic, is of mi- 
nor importance. Note however that the classical variance 
a determines the height and the correlation field ipc- 
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FIG. ^. Contributions to the susceptibility for kpL = 60 
and £/L — 2, with (normalized) semiclassical estimates for 
(x'^iO)) (dashed-dotted line) and (x^(0)) (short dashed) for 
a generic chaotic geometry and for (x'*(0)) (solid line) and 
(x^(0)) (long dashed) for square geometry, xo ~ 1xl| for the 
square and \xl\L^ /96cy for the chaotic geometry. 

The susceptibility (x'*(0)) a-s a function of for 
£/L = 2 is shown as the dashed-dotted line iiL.Fig. ^ 
and a corresponding approximation for (x^(0))cJ for fi- 
nite Lrh, 



(xHo)) 



96crLrf 



\Xl\ 



(25) 



is shown as the short dashed line in Fig. |^. The two con- 
tributions ( p^ , p5| ) add up to an £- independent magnetic 
response 



(X(0)) 96aL^ 



\Xl\ 



L4 



(26) 



which is the same as that for clean chaotic systems, given 
the assumption of an ^-independent variance. 

For comparison, results for the square geometry as a 
function of for £/L = 2 are also presented in Fig. |^. 
The solid line shows our numerical results for {x'^{0)) and 
the long dashed line shows (x^(0)) (see RehEl for more 
details). The dependence also shows qualitative sim- 
ilarity between the disorder-induced susceptibility in the 
square and the chaotic geometry. Naturally, the square 
and the chaotic billiard show quite different behaviour 
for the size-induced susceptibility. 

For the square a crossover from domination by 
disorder-induced correlations to domination by size- 
induced correlations occurs for ~ which contrasts 
with chaotic geometries where the crossover occurs for 
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V. CONCLUDING REMARKS 

We considered disorder-induced spectral correlations of 
mesoscopic quantum systems in the non-diffusive regime. 
We combined a diagrammatic treatment of the disorder 
with a semiclassical approach for the Green functions 
involved. This leads to the representation of the dia- 
grams involved in terms of (pairs of) classical paths in 
real space which undergo both (multiple) scattering at 
the impurities and reflections at the system boundaries. 
This hybrid-type approach enables us to perform disorder 
and energy averages and to study the complete crossover 
from the diffusive to the clean limit for arbitrary values of 
the elastic mean free path t and (smaller than vptH)- 
In particular it is applicable to ballistic micro-structures 
of, in principle, arbitrary geometry. The approach, pre- 
sented here for billiard systems, applies also to systems 
with potentials whose effect is then incorporated in the 
semiclassical Green functions. 

We focussed on the effect of disorder-induced corre- 
lations on the averaged magnetic susceptibility which is 
closely related to the spectral number variance and which 
amounts to evaluate Cooper type contributions to the 
correlator. Corresponding diffuson type diagrams can be 
computed equivalently. As an example of a system with 
integrable clean counterpart, we treated in detail the (ex- 
perimentally relevant) square billiard. We showed that 
there are two distinct regimes of behaviour depending on 
the relative magnitudes of £ and L^. Certain statisti- 
cal assumptions on the trajectories in ballistic systems 
with a generic chaotic geometry enabled us to derive ap- 
proximate analytical expressions for the disorder induced 
magnetic response. It turns out that the magnetic field 



dependence and L^-dependence of the susceptibility, at 
least qualitatively, resembles that of the integrable case. 
This implies that the actual geometry of the clean system 
is of minor importance for the disorder induced orbital 
magnetism even in the ballistic regime. 
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APPENDIX A: DIAGONALISATION OF C 
OPERATORS FOR ZERO FIELD 

In this appendix we show that within the semiclassical 
approach the analytical diagonalisation of the Q opera- 
tors, Eq. (^, for a square geometry at zero field is pos- 
sible and straightforward. We employ the extended zone 
scheme which is constructed by reflecting the original bil- 
liard with respect to its boundaries. Any given trajectory 
connecting ri = {xi,yi) and r2 = {x2,y2) involving an 
arbitrary number of specular reflections at the bound- 
aries is tranformed into a straight line path connecting ri 
and y\ where y\ = (±a;2 -t- in^L, ±j/2 + '2nyL) is an image 
of the position r2 = {x2, 2/2), and Uy are integers. As a 
result a semiclassical propagator, C'^'^H^^i j Z/i! ^2,1/2), may 
be written as a sum of propagators, Q'"^xi,yi;x2,yl), 
along straight-line trajectories, 



J 



C^'^\xi,yi;x2,y2) = C^^\xi,yi\±X2 + 2nxL,±y2 + 2nyL). 



(1) 



Since C,'^'-^^ is dependent only on the distances X = {x2 — xi) and Y = (1/2 — yi) along a straight line it is possible to 
introduce a Fourier transform X{qx, %) and, using the Poisson summation rule, we may write 



C,^^\xi,yi;x2,y2) = ^ J '^^2tt)'^ 4: X{qx, Qy) cosiq^xi) cos{qyyi)d f — -ma;^ S (^^-niy | e 



(2) 



Performing the integrals gives quantised momentum values = mxTr/L, qy — myn/L and 



C^'^Hxi,yi;x2,y2) = 



OC 00 



A(0,0)+4^ E 



I cos(- 



niT — l m., — 1 



^ -)cos( ^ )cos(- ^ 



-) cos(- 



+ 2Y: A(^,0)cos( 



m^ — l 



^rrixTTXi ,mxTrx2 
) cos( 



) + 2 ^ A(0,^)cos(^)cos(!^) 



(3) 
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Although (■('-^^ {xi,yi; X2,y2) is not translationally invari- 
ant, it does have a periodicity of 2L. The diagrams Sn 
may be written in terms of the Fourier transform \{qx, Qy) 
by inserting the above expression into Eq. (^). Spatial 
integrals are performed using the following relation, 



^ 



and the result is 



(4) 



(5) 



m-c—O rrtoj—O 



In order to evaluate X{qx,qy) we Fourier transform the 12 



expression for the semiclassical operator Q , Eq. (15), 
performing integrals over X and Y. By changing from the 
cartesian coordinates {X,Y) to polar coordinates {r,9), 
the radial integration is done giving 



K<lx,qy) = 
2^ 



(6) 



dO 



/Q (1 + 7T — zcjr) + i£ [qx cos + < 
and we find 

1 



, sm( 



Klx,qy) 



{\ + -IT - iujrf + {tqf 



(7) 



where q^ 



+ q'y- This result is iri agreement with 
the approximate diagonalisation of Ref.u. Together with 
Eq. it gives the spectral correlation function in the 
"ballistic" regime. 
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